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THE Za INDEX OF DISORDERED TOPOLOGICAL 
INSULATORS WITH TIME REVERSAL SYMMETRY 

HOSHO KATSURA AND TOHRU KOMA 


Abstract. We study disordered topological insulators with time re¬ 
versal symmetry. Relying on the noncommutative index theorem which 
relates the Chern number to the projection onto the Fermi sea and 
the magnetic flux operator, we give a precise definition of the Z 2 index 
which is a noncommutative analogue of the Atiyah-Singer Z 2 index. We 
prove that the noncommutative Z 2 index is robust against any time- 
reversal symmetric perturbation including disorder potentials as long as 
the spectral gap at the Fermi level does not close. 


1. Introduction 

In a seminal paper, Haldane m proposed a tight-binding model of spin¬ 
less fermions on a honeycomb lattice. Interestingly, although the net mag¬ 
netic flux through the unit cell of the lattice is vanishing, the model pos¬ 
sesses a topological phase in which the energy bands carry nontrivial Chern 
numbers [33 EH- Kane and Mele [23] generalized the Haldane model to a 
spin-I/2 fermion model with spin-orbit coupling, so that the model is time- 
reversal invariant. Because of the time-reversal symmetry, the sum of Chern 
numbers for occupied bands is always vanishing. Therefore, at first glance, 
it appears that the model is topologically trivial. Kane and Mele, however, 
showed that the model exhibits a Z 2 topological order which is characterized 
by the Z 2 topological invariant, instead of the Chern number. Further, their 
study indicates that the Z 2 topological invariant they constructed is related 
to the Z 2 index which was introduced by Atiyah and Singer mm. in order 
to characterize topological properties of a manifold. The argument of Kane 
and Mele relies on the properties of the Bloch wavefunctions, and hence 
requires translational invariance. Therefore, one of the important issues is 
to find a mathematically rigorous and numerically accessible representation 
of the Kane-Mele Z 2 index which is robust against perturbations such as 
disorder potentials^ 

We should also remark that Schulz-Baldes [33 Illlg] pointed out that 
the Z 2 index by Kane and Mele is equal to the Z 2 index of the Fredholm 
operator which is written in terms of the projection onto the Fermi sea 
and the magnetic flux operator. As is well known, in the case of the stan¬ 
dard quantum Hall system, the usual index of the corresponding Fredholm 


Date. March 3, 2016. 

^ Following the Niu-Thouless-Wu argument [28] , Kane and Mele proposed a generaliza¬ 
tion of the Z 2 index to systems with disorder and interactions. However, the robustness 
and the Z 2 quantization of such an index have not yet been proved rigorously. 

1 



2 


H. KATSURA AND T. KOMA 


operator is equal to the Chern number. This relation is nothing but the non- 
commutative index theorem [asKi] which is a noncommutative analogue 
of the Atiyah-Singer index theorem. In general, such a Fredholm operator 
is noncompact and its spectrum contains the essential spectrum. Thus, it 
seems hard to prove the continuity of the eigenvalues of the noncompact 
Fredholm operator against generic perturbations such as disorder, although 
the continuity leads to the robustness of the Z2 index. 

In the present paper, instead of the Fredholm operator [Ml He] , we use 
a pair of projections which was introduced by mis]. More precisely, we 
give an operator-theoretic definition of the Z2 index by relying on the non¬ 
commutative index theorem mum which relates the Chern number to the 
index of the pair of the projections. One of the projections is the projection 
onto the Fermi sea and the other is its unitary deformation which physically 
implies the insertion of a magnetic flux into the system mi- The advantage 
of our approach is that the difference between the two projections becomes 
a compact operator. Therefore, all of the spectrum are discrete with fi¬ 
nite multiplicity except for zero. This makes the treatment of perturbations 
much easier than that for noncompact Fredholm operators. 

In conventional differential geometry mm, the Z2 index is defined by the 
analytical index of an operator such as an elliptic pseudo-differential opera¬ 
tor on a manifold. A special property of the multiplicities in the spectrum 
of the operator yields the robustness of the Z2 index [ 2 ]. On the other hand, 
in our approach based on noncommutative geometry |H], the corresponding 
Z2 index is defined by the analytical index of the above compact operator, 
say A, consisting of the projection onto the Fermi sea and the magnetic flux 
operator. Similarly to the Kramers degeneracy in the energy spectrum, the 
time-reversal symmetry leads to a special property of the degeneracy in the 
spectrum of A: Every eigenvalue A of A shows an even degeneracy when 
0 < |A| < 1 . This property yields the robustness of the constructed Z2 index 
as in the case of the Atiyah-Singer Z2 index. 

The present paper is organized as follows: In the next section, we de¬ 
scribe the models we consider. In Sec. 3 , we state and prove our main 
result (Theorem [ 31 ) about the noncommutative Z2 index. As an example, in 
Sec. m we demonstrate that the Z2 index for the Kane-Mele model can be 
computed from the Chern number for the corresponding spinless Haldane 
model. Appendix [Mis devoted to a short review of noncommutative geom¬ 
etry, in which the relation between the analytical index of an operator and 
the Chern number is given. In Appendix [Bj we present a calculation of the 
Chern number of the projection onto the Fermi sea for the Haldane model 
in the present setting. In Appendix [Cl we demonstrate that the present Z2 
index is equal to that by Kane and Mele for the translationally invariant 
system. 


2. Models 

Consider a tight-binding model of spin- 1/2 fermions on the square lattice 
Z^. We denote by f and j, the spin-up and the spin-down indices, respec¬ 
tively. We assume that there are r orbitals at each lattice site. For each 
site n G Z^, we denote by (^^(n, /x) G C the amplitude of a wavefunction 


THE Z 2 INDEX OF DISORDERED TOPOLOGICAL INSULATORS 


3 


Lp G C'" (8* C^), where a G {t, i} and /x = 1, 2,r denote the spin and 

the orbital indices, respectively. The action of the Hamiltonian H on the 
wavefunction ip is written as 

r 

( 2 . 1 ) {H(p)c,in,p) = 

mGZ2 u=l /3G{t,l} 

for n G Z2, p, = 1,2,... ,r, and a G {t)'!-}) where are the hoping 

integrals satisfying the hermitian condition, 

_ j.P,a 

Here, stands for the complex conjugate. In the following, we assume 
that the range of the hoping integrals is finite. The two-component spinor 
wavefunctions for the p-ih. orbital at the site n G is written in the form 

Let US introduce the time reversal transformation 0 for wavefunctions as 

(2.2) tp® := Otp = for (^g£2(z2,C’’®C2), 

where t/® is a unitary transformation which can be written as a product 
of local unitary transformations with a compact support, and has a finite 
period, i.e., invariant under a finite translation on the present lattice Z^. As 
usual, the complex conjugates for wavefunctions are defined by 

{Tp)ain,p) = ipaiT^,p) for n G Z^, /x = 1, 2,..., r, and aG{t,i}- 

In what follows, we assume that the time-reversal transformation 0 satisfies 

0V = -^ for 

and say that 0 is an odd time-reversal transformation. 

Let a be an operator on the Hilbert space £^(Z^, C’’(8>C^). We say that the 
operator a is odd time-reversal symmetric if the following condition holds: 

(2.3) 0(a(p) = a(p® for any (p G £^(Z^,C^ (g) C^). 

3. The Z 2 Index 

We assume that the Fermi level lies in a spectral gap of the present Hamil¬ 
tonian H. We write Pp for the spectral projection on energies below the 
Fermi energy Ep. We introduce a unitary transformation Ua for wavefunc¬ 
tions ip G (g) C^) as 

{Ua‘p){'n, p) = Ua{'n)ip{n, p) for n G Z^, and ^=1,2, ...,r, 

where 

(3.1) C/.(u) := 

|xxi -I- IU 2 - (oi -I- ia 2 )\ 

for u = (xxi,xx 2 ) G Z^ and a = ( 01 , 02 ) G (Z^)*. Here, (Z^)* denotes the dual 
lattice of Z^, i.e., (Z^)* := 7? — (1/2,1/2). 

We write 

(3.2) A = Pp- UaPpU* 
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for short. Clearly, this operator A is the difference between two projections. 
As in mum, the operator is trace class0 and the relative index can be 
defined as 

(3.3) Ind(PF, UaPpU*) := dim ker (A — 1) — dim ker (A + 1) = Tr A^, 

where dim ker O stands for the dimension of the kernel of an operator O. 
This implies that the quantity Tr A^ is quantized to an integer. Actually, 
this is nothing but the Chern number. See Appendix!^ for details. 

Consider a Hamiltonian H of the form m which is odd time-reversal 
symmetric. We will show that the Z 2 index for such a Hamiltonian H can 
be defined by 

(3.4) Ind 2 (PF, U^PpU*) := dim ker (Pf - U^PpU* - 1) modulo 2. 

More precisely, we will show that the Z 2 index constructed above is robust 
against any perturbation which has the same odd time-reversal symmetry as 
the unperturbed Hamiltonian whenever the Fermi level lies in the spectral 
gap of the total Hamiltonian. 

We write 

(3.5) b = 1-Pf- U^PfU* 

for short. Then, as in [5], the following two relations hold: 

(3.6) AB + BA = 0 and + B"^ = 1. 

Here, the operator A is given by ra . 

We note that the spectrum of the operator A of (13.2p is discrete with 
finite multiplicity except for zero because A is self-adjoint and A^ is trace 
class. 

Following [5], we use the standard supersymmetry argument to prove 
that many of the eigenvectors of the operator A come in pairs related by 
the operator B. Let y? be an eigenvector of A with eigenvalue A G (0,1), 
i.e., Aip = A(/?, with A G (0,1). Then, the anticommutation relation of (13.6p 
yields 

ABif = —BAip = —XB(p. 

Further, the second relation of (13.6p yields 

hV = = (1-a^)v?. 

These results imply that B is an invertible map from an eigenvector with 
eigenvalue A G (0,1) to that with —A. Thus, there is a one-to-one corre¬ 
spondence between the states (p and B(p, and their eigenvalues come in pairs 
±A, provided that 0 < |A| < 1. 

Next, we show that the time-reversal transformation 0 plays a similar 
role as the operator B in the above argument. Let (p be an eigenvector of A 
with a strictly positive eigenvalue, i.e., Acp = Xip, with A > 0. From (|3.2I) . 
one has 

(3.7) 0(Pf - U^PFU*)ip = Xip^. 


^For trace class operators, see, e.g., Chap. VI.6 of the book m- 
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Since the Hamiltonian H is time-reversal symmetric, the spectral projection 
Pp is time-reversal symmetric, too. We assume that the unitary operator 
[7® in the time-reversal transformation 0 of (12.2p satisfies 

(3.8) U^U^iU^y = Pa- 
Under these assumptions, Eq. (1321) can be rewritten as 

(Pf - P^PpPa)^® = A<^®. 

Further, the left-hand side of the above equation is written as 
U*{U^PfK - -PF)Pa<^® = -U*AU^ip^. 

Consequently, we have 

H(Pa(/P®) = -APa(/^®. 

This implies that PaV?® is an eigenvector of A with eigenvalue —A. Further¬ 
more, the transformation is invertible. 

Lemma 1. Let cp be an eigenvector of A with eigenvalue A satisfying 0 < 
|A| < 1. Then, the maps, B and Pa©, are commuting with each other when 
acting on the vector ip, i.e., 

(3.9) U^iBipf = P(Pa¥^®). 

Proof. From the definition (j3.5p of the operator B and the assumption (|3.8p , 
one has 

e{Bip) = (1 - Pp - u*PFU^)ip^ 

= U*{1 - PaPpP: - PF)Pa<^® = KBU^ip^. 

Therefore, we have (|3.9n . □ 

Clearly, two vectors cp and Ua{Bip)® are eigenvectors of the operator A 
with the same eigenvalue A. If these two vectors are linearly independent 
of each other, then the corresponding sector which is spanned by the two 
vectors is a two-dimensional space. 

Lemma 2. Let cp be an eigenvector of A with eigenvalue A satisfying 0 < 
|A| < 1. Then, 

(3.10) {ip,U^{Bipf) = 0 

Proof. For a vector if G C'’ ® C^), one has 

(3.11) {ey,Q(p) = {-ip,ip) = {(p,y) 

from the definition of the time-reversal transformation 0. We set 

y = uyBipf. 

Substituting this into the right-hand side of the above equation (|3.1ip . we 
have 

(3.12) {ey,eip) = {ip,uyBipf). 

On the other hand, the above vector if can be written as 

y = Q{KBp>), 
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in the same way as in the proof of Lemma [TJ Substituting this into the 
left-hand side of Eq. (j3.12l) , we obtain 

(3.13) - Qif) = ( 99 , 

where we have used the fact that = —(p for any vector cp. Then, we find 
that the inner product in the left-hand side can be written as 

{KBip,eip) = {<f,BiU,ip^)) = {ip,U,iB<ff), 

where we have used the relation (j3.9|) in Lemma[T]to get the second equality. 
Because of the minus sign in the left-hand side of (|3.13p . we arrive at the 
desired result (j3.10|) . □ 


These two lemmas imply that the multiplicity of the eigenvalue A of the 
operator A must be even when A satisfies 0 < |A| < 1. Therefore, if the 
eigenvalues of A change continuously under a continuous variation of the 
parameters of the Hamiltonian H, then the parity of the multiplicity of the 
eigenvalue A = 1 must be invariant under the deformation of the Hamilton¬ 
ian. Thus, we have two possibilities: (i) An even number of eigenvectors 
of A are lifted from the sector spanned by the eigenvectors of A with the 
eigenvalue A = 1; (ii) an even number of eigenvectors of A with eigenvalue 
A 7 ^ 1 become degenerate with the eigenvectors of A with the eigenvalue 
A = 1. In both cases, it is enough to prove the continuity of the eigenvalues 
of the operator A under deformation of the Hamiltonian [321E2], in order 
to establish our main result. 

Let us consider a perturbation 5H for the odd time-reversal symmetric 
Hamiltonian H. We assume that the perturbation 5H is odd time-reversal 
symmetric with respect to the same time-reversal transformation 0 as the 
unperturbed Hamiltonian. We also assume that the range of the hopping 
integrals in the perturbed Hamiltonian 6H is finite. Therefore, the norm 
||(5iL|| of the Hamiltonian 6H is finite. Clearly, the total Hamiltonian H' = 
H+6H is odd time-reversal symmetric, too. We assume that the Fermi level 
of the Hamiltonian H' still lies in the spectral gap of H'. Using a contour 
integral, the projection onto the Fermi sea for the Hamiltonian H' is written 


(3.14) 


Pf 


1 

27ri 


dz 


z-H'' 


where the contour encloses all the spectrum of H' below the Fermi level. 
Similarly to the operator A of (13.2p . we introduce 


A'= P'- U^PLU: 


Then, the difference between A and A' is written as 


A' -A = {P^- Pf) - C/a(FF - P^Wl 

In order to evaluate this, it is enough to estimate Pp — Pp. Using the contour 
integral, we have 


Pf-Pf = 


2TTi 


dz 


1 


1 




z-H' 
1 


SH 


z-H 

1 


z-H' z-H 
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Because of the assumption of the spectral gap, the norm of the right-hand 
side can be bounded by the norm of the perturbation 6H with some positive 
constant. Therefore, the difference A' — A is bounded by the norm, too. In 
consequence, the operator A is continuous with respect to the norm of the 
perturbation 5H. Using the min-max principle0 we obtain the desired result 
that the nonzero eigenvalue A of the operator A is continuous with respect 
to the norm of the perturbation 5H. 

To summarize, we obtain: 

Theorem 3. Suppose that the Hamiltonian H is odd time-reversal symmet¬ 
ric with respect to the transformation 0 of / fH.gj) whose unitary operator U® 
satisfies the condition id. (S|) for the unitary operator Ua of the operator A of 
19. Sf) . We assume that the Fermi level lies in a spectral gap of the Hamilton¬ 
ian H. Then, the Z2 index of \3.4^ is continuous with respect to the norm 
of any perturbation if the perturbed Hamiltonian preserves the same time- 
reversal symmetry as the unperturbed Hamiltonian. In other words, as long 
as the spectral gap does not close, the Z2 index is robust against any pertur¬ 
bation which preserves the same time-reversal symmetry as the unperturbed 
Hamiltonian. 


Remark. (i) When the Fermi level lies in the localization regime, the same 
statement is still valid. But, we need to prove the localization of the particles 
separately. More precisely, we need a decay estimate for the resolvent {z — 
H)~^, z G C. See, e.g., [32l [22] for the homotopy argument in such a 
case. See also [22] which discusses the conditions under which both the 
quantization and homotopy invariance of the Chern number in four or higher 
dimensions hold in the presence of strong disorder. 

(ii) Since the Z2 index is given by the dimension of the kernel of the operator 
as in dsai), one might think that the Z2 index may be written in terms of 
an integral of some connection, similarly to the Chern number. In fact, the 
Z2 index which was defined by Kane and Mele can be written in terms of an 
integral of the same connection as that of the Chern number over one-half 
of the Brillouin zone for translationally invariant systems [1211261 ITT] . (See 
also related articles [laiMiiiaisi].) The explicit computation of the Z2 
index with the use of the integral formula is given in Appendix [Cj 

(iii) As mentioned in Introduction, Schulz-Baldes [Ml [IS] defined Z2-indices 
for general odd symmetric Fredholm operators and treated the Z2 index for 
the Kane-Mele model as an example. (See also |141 (9] for related articles.) 
But his approach is different from ours. In fact, he defined the Z2 index by 
the parity of dim ker T for the Fredholm operator T := PpUaPp + (1 ~ .Pf)- 
From Eq. (3.2a) in the proof of Proposition 3.1 in |5], the dimension of the 
kernel of T coincides with ours as dim ker T = dim ker (A — 1) for the 
operator A of ()3.2p . On the one hand, the operator T is noncompact and 
contains the essential spectrum. On the other hand, the operator A has 
only the discrete spectrum with finite multiplicity except for zero, because 
A is compact which follows from the fact that is trace class. Therefore, 
a homotopy argument for A is much easier to handle than that for T. It 

^For the min-max principle, see, e.g.. Chap. XIII.1 of the book of m- 
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should be noted that the homotopy argument is indispensable for defining 
the Z 2 index for disordered systems. 

(iv) Hastings and boring usmaES] proposed an alternative noncommuta- 
tive approach to the Z 2 index of topological insulators. 


4. The Z2 Index of the Kane-Mele Model 

As a nontrivial example, we demonstrate in this section that the Z 2 index 
for the Kane-Mele model can be computed from the Chern number for the 
corresponding Haldane model with a specific flux pattern and vanishing 
staggered sublattice potential. Both Hamiltonians of the Kane-Mele and 
the Haldane models have a single spectral gap except at critical points. In 
what follows, we assume that the Fermi level lies in the spectral gap. We 
first note that the condition (j3.8p holds for the Kane-Mele model. Because 
of the continuity of the Z 2 index, it is sufficient to treat the simplest case 
without either disorder or the Rashba term in the Hamiltonian. Namely, 
when continuously switching on those terms, the Z 2 index does not change 
as long as the spectral gap does not close. As is well known, the Hamiltonian 
of the Kane-Mele model is decoupled into two independent Haldane models 
HZ! when the Rashba term in the Hamiltonian vanishes. Although each of 
the Haldane model shows the nonvanishing Chern number, the total Chern 
number for the decoupled Haldane models is vanishing because of the time- 
reversal symmetry. This can be seen as follows. First note that the one-to- 
one mapping sends an eigenvector of the operator A with the eigenvalue 
A = 1 to that of A with A = — 1. Then, since the difference between the 
multiplicities of A = 1 and A = — 1 gives the total Chern number, it follows 
that the total Chern number is vanishing. 

By contrast, the Z 2 index is given by the parity of the multiplicity of 
A = 1. Therefore, from these observations, we reach a conclusion; When the 
Rashba coupling is weak so that the spectral gap does not vanish, the Z 2 
index for the Kane-Mele model is equal to the parity of the Chern number 
for either one of the decoupled Haldane models. Namely, we have 

(4.1) lnd 2 {PF,UaPFU*) = Ich modulo 2, 

where we have written for the Chern number for the single Haldane 
model. In Appendix iBl we show that Ich = 1 for the single Haldane model 
in the present setting. Therefore, the Z 2 index for the Kane-Mele model is 
equal to 1. 


Appendix A. The Index Theorem and the Chern Number 

Using the method of noncommutative geometry isiEiiaii], one can show 
that the quantity, Tr A^, in the right-hand side of the index formula (13.311 
is equal to the Chern number. In this appendix, we do not require that the 
Hamiltonian H is odd time-reversal symmetric. 

By using a function 0a (u) which is the angle of sight from a to u, the 
unitary operator f7a(u) of (13. ip can be written as t/a(u) = exp[i0a(u)]. We 
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introduce an orthonormal complete system of functions, 

tv. — n = u, /3 = a and u = /r; 

“ ^ ’ | 0 , otherwise. 

Namely, the function Xu^ has the nonvanishing value 1 only for the spin a. 
and the /r-th orbit at the site u. One has the matrix elements, 

where 

Tu,v = 1 - exp[i6'a(u) - i6'a(v)] 

and 

(A.i) 

Using the complete basis, the index of (|3.3p is written 
Ind(PF, = Tr (Pf - 

— ^ ^ ^U,V^V,W^W,uSu,V,W,U; 


where 

(A.2) 




Ou,v,w,u — -'u.V ‘'v.W ‘'WjU 

As shown in Proposition 3.8 in [1], the index Ind(PF; UaPpU*) is independent 
of the positioning a of the flux tube. Relying on this fact, the index is written 

IndpPp, C/aPFUa) — 1.^1 ^ ^ ^ ^ 'Fu,vTv,wTw,uSu,v,w,U) 

aeA* u,v,weZ2 

where A* is a finite rectangular box which is a subset of the dual square 
lattice (Z^)*. Following Elgart, Graf and Schenker [10] . one has 

Lemma 4. There exists a sequence of finite rectangular boxes A C such 
that 

Ind(PF ) UaPFUa) — ^-^2 ^ ^ ^ ^ ^ ^ Tu,vTv,wTw,uSu,v,w,U' 


uGA v,wSZ 2 aS 


Further, in order to compute the right-hand side of the index, we apply 
the Connes’ area formula 13 S 1 II], 

Tu,vTv,wTw,u = 27ri(v - w) X (w - u) 


aG 

= 27ri[(u2 - W2){wi - Ul) - (ui - Wl){w2 - U2)\. 
In consequence, the Index is written 

Ind(PF, UaPpU*) = lim 777 Y] (v - w) x (wr - u)Su v,w,u- 

Atz2 A ^ ^ 

' ' ugAv,wGZ2 

Note that 

(U 2 -u; 2 )(rci -Ul)(xr,^FX^’")(x^’^^FX;^^)(x;^^^’FXr) 

= (xr,^FXe’")(xe’^ [X2,Pf]xI^){xI^, [Ai,PF]xr) 
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and 

(ni -ini)(in2 -n2)(xr,^FX^’")(xe’^^’FX;^«)(x;^^^’FXr) 

= (xr,^Fxe’")(xe’^ [XuPf]xI^){xI^, [^2,PF]xr), 

where the operators Xj are given by (Xjy?)(n) = nj(p(n), j = 1,2, for 
(/?(n) G C’’ ® C^. Combining these, (|A.ip and (IA.2I1 . we obtain 

lnd{PF,U^PFU*) = -2m Im ^Tr xa^^f[[^i, ^f], [X2,Pf]]xa, 

|A| 

where XA is the characteristic function of the finite lattice A, i.e., 

1, n G A; 

0 , otherwise. 

Further, the index can be expressed as miiniiss] 

Ind(PF, U^PfK) = -27rfTr PF[[l?l,a, ^f], [^? 2 ,a,PF]] 
in terms of the step functions, 

1, rij — aj > 0; 

0, rij — aj < 0, 

for j = 1, 2. Using a contour integral as in (j3.14p . we have 


XA(n) = 




[«i...ai = A 

1 

27ri 

i 

2'Ki 

where we have introduced 


1 


1 




dz- 


1 


dz 


z-H 
1 




'z-H 


J. 


z-H ^'‘^z-H’ 


the local current operators Jj^a = i[H,'dj^a] 
for j = 1, 2. Substituting this into the expression of the index, one has 


(A.3) Ind(PF,C/a/^FC/a*) 
1 


Xfdzf dz' Tr 


-H 


-(lo 2 ). 


Consider the system on the finite lattices A. When there exists a uniform 
spectral gap above the Fermi level with respect to the size of the finite 
lattices A, the index is written as 


(A.4) Ind(PF,i7a^’FC/:) 


= 27rz lim 
AtZ 2 


A A 

i:Ei<CEY j-Ej'>Ep 


1 

(Ei - Ej)^ 


J 2 ,a^i) -( 1 ^ 2 ) 


in terms of the eigenvectors of the finite-volume Hamiltonian H with the 
energy eigenvalue Ei. Here, Ef is the Fermi energy. The right-hand side 
expresses the Hall conductance which is quantized to an integer because the 
analytical index Ind(PF; U&PFUt) of operator in the left-hand side takes 
an integer value by definition. 
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Appendix B. The Chern Number of the Haldane Model 


In this Appendix, we present a detailed exposition of the calculation of 
the Chern number for the Haldane model. The original Haldane model is 
defined on a honeycomb lattice. For technical reasons, however, we need to 
define it on the square lattice T? in the present setting]^ For this purpose, 
we introduce two orbitals labeled by a and b at each lattice site n G Z^. We 
write 


<^(n) = 


G 


for the wavefunction if G at the site n = (ni,n 2 ). In order to 

prove the statement in Sec. 01 it suffices to consider the Haldane model with 
a pure-imaginary hopping term between second-neighbor sites and the van¬ 
ishing staggered sublattice potential. In this case, the Schrodinger equation 
reads 


t[ip^{ni - 1, n2) + 112 ) n2 1)] 

-F - 1,71,2 - 1) - ip°'{ni,n2 - 1) + -h 1, 71 - 2 ) - + I,n2 + 1) 

-F (^“(ni, 712 + 1) - - 1, 712 )] = Eip'^{ni,n2) 

and 


712) + V 9 “(ni, 712 - 1 ) + + 1 ,712)] 

-I- it'[ip’^{ni - 1, 712 ) - (p^{ni — 1,712 - 1) -I- ip^{ni,n2 - 1) - f^{ni + 1, 712 ) 
-I- ip’^{ni -I- 1,712 + 1) - </^^(iii, 112 + 1)] = Fi(/3^(?ii, 712 ), 


where the hopping integrals, t and t', are real constants, and E is the energy 
eigenvalue. In the following, we will treat only the case with t' > 0. 

In order to use the formula ()A.4p for computing the Chern number, we 
consider a finite Li x L 2 rectangular box with the periodic boundary con¬ 
ditions in 7?. We set 


V9(n) 




1 

VL 1 L 2 


exp[zk • n] 




where k = {ki, ^ 2 ) is the wave number vector, and and are functions 
of k. Substituting this into the above Schrodinger equation, one has 



where 


with A(k) 


H(k) 


2t'[sin ki + sin k 2 


/-A(k) F(k)\ 

V r(k) A(k); 

sin(A:i -|- /C 2 )] and 


(B.l) 


r(k) = t(l + e*^i 


The energy eigenvalues are then given hy E = Fi±(k) = ±£(k) with £(k) = 
VA(k)2 + |F(k)|2. 


^ Note that an arbitrary two-dimensional tight-binding model can be mapped onto the 
model on 1? with suitably chosen hopping integrals. 
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Let US consider the case of r(k) = 0. In the Brillouin zone, (—7r,7r] x 
(—TT, vr], there are only two wave-number vectors that satisfy this condition. 
They are given by {ki,k 2 ) = ±(27r/3, 27r/3). For these two points, one has 
A(k) = ±3\/3t^ for the positive and negative wave numbers, respectively, 
because we have assumed t' > 0. Thus, there exists a nonvanishing spectral 
gap between the upper and lower energy bands for t 7^ 0 and t' / 0. 

The eigenvector of the lower band is given by m 

^ 1 /L(k) + A(k)\ 

V^2£(k)[£(k) + A(k)] V -r(k) ;■ 

But, the normalization factor of this vector becomes infinite at {ki,k 2 ) = 
— (27r/3, 27r/3) because A(k) = —£(k) at this point. An alternative expres¬ 
sion of the eigenvector is given by m 

o fki = 1 ( -m 

v'2£(k)[£(k) - A(k)] V£(k) - A(k)^ 

In this case, the normalization factor of this vector becomes infinite at 
{ki,k 2 ) = {2tt/3,2tt/3). The relation between these two vectors is given by 
the gauge transformation, /-(k) = g'_ (k) , for {ki, /C2) / =t(2vr/3, 27r/3), 

where the angle function 17 (k) is given by 

^*r,(k) _ r(k) 


(B.2) 

We choose the eigenvector as 


|r(k)r 




and 


^'°"(k) = 


<^'-(k) 
<^'-(k) 


=/_(k) for (/ci, A:2) G := (— TT, tt] X [0, vr] 


= ^-(k) for (fci, A:2) G := (— TT, tt] X [—TT, 0]. 


Using these eigenvectors and the translational invariance, the Chern number 
Ich which is given by the right-hand side of (|A.4D is written as 


(B.3) 

with 

(B.4) iL"P = 




low 


/ ^ 


dkidk2 


and 

(B.5) 


^low _ 




dkidk2 


(^"p(k), Ji(k)^:t"(k))(^:"(k), J2(k)^:ip(k)) 

[i?+(k)-E_(k)]2 

- (1 GA 2) 

((^'°»(k), Ji(k)(^lp-(k))((^'P«(k), J2(k)(^l°-(k)) 


[F;+(k) - i?_(k)]^ 


-(1 ga2 ), 


where (^!|P(k) and <^!f!"(k) are the corresponding eigenvectors for the upper 
band, and the current operators are given by 
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By differentiating 

(<^:iP(k),/7(k)^7(k)) = o, 

one has 

+ (<^!?(k), ■^H{k)ip'^{k)) 

+ ((^!?(k),iJ(k)^(^7(k)) = 0. 

From this, the definition of the current Jj{k) and the identity, 

(B.6) (^(,5“P(k),(,57(k)) + ((^"’’(k), ■^<f'^{k)) = 0, 

one obtains 

((^"’’(k),7(k)(^7(k)) = [F;-(k) -F;+(k)](^(^"P(k),(^7(k)). 

In the same way, 

= [E-{k) -F;+(k)]((^7(k),^(^"P(k)). 

Substituting these into the integral K'^'^ of (IB.4D and using the identities 
(IB.6I) and 


4-V®(k)) = 0, 


'dki 


dki 


one has 


(B.7) 


= ^^pdA;idA;2(^(^!iP(k),(^7(k))((^7(k),^(^"P(k)) - (1 o 2) 


/ 

Ju^ 

I 

Ju^ 

/ 71 

-1 


dkidk2 


dkidk2 


dki 


^W-(k) A^-(k))-A(^-(k).A^^p(k)) 


In the same way, 

(B.8) = -I" dfci(7r7i,0),^7-"'(^i,0)) 

+ j dki{(p^°^{ki,-7T),-^ip^^'"{ki,-7T)). 


On the other hand, one has 


uki 


drj{k) 

dki 


7T(k) + e*^(*^)T^7T(k). 

uki 
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from the gauge transformation, Combining this 

with ()B.7p and (IB.8D . one obtains 

/ TT O /* 7^ 

dki—7]{ki,0)-J dki—r]{ki,Tr) =-27ri, 

where we have used (jB.ip and ()B.2D . Substituting this into equation (IB.3D . 
the value of the Chern number is calculated as Ich = 1 for the present 
Haldane model. 


Appendix C. Relation to the Z 2 index by Kane and Mele 


In order to show that the Z2 index by Kane and Mele is equal to the 
present index (I3.4|) . we recall the integral formula |12j for the Z2 index. We 
use the formula given by Eq. (9) in [TT], i.e., 

(C.l) D:=^ I dk-A- [ dkidk2 3^ , 

27! ya(EBZ) JeBZ J 

where A and 3” are, respectively, the Berry connection and the field strength, 
and the “effective Brillouin zone”, EBZ, stands for one-half of the Brillouin 
zone with a time-reversal invariant frame. Because of the homotopy argu¬ 
ment, it suffices to consider the case where the Kane-Mele model is decoupled 
into two independent Haldane models. 

Consider first one of the Haldane models whose lower band carries the 
Chern number -|-1. Hollowing the argument in Appendix [Bl we set EBZ = 
1X“P. Further, we choose 


Vk5'-(k)) 

as the connection A in the line integral in (jC.lh . For the connection of the 
field strength 3", we choose 


yi;(k) = -i(/_(k),Vk/-(k)). 

Here, we stress that the wavefunction /_(k) has no singularity in U“p, while 
the wavefunction (7-(k) has the single singular point in IX^p. Then, the 
contribution to the index D is written as 


D+ := 


1 


® dk ■ Ag 

JdU'^p 


® dk ■ Af 

JdU'^P 


1 


/ 

JdU'^P 


dk ■ Vki/(k) = 1 


in the same way as in Appendix [Bj 

Next consider the other Haldane model which is the time reverse of the 
above one and therefore whose lower band carries the Chern number — 1. 
The connection yi(k) in the line integral can be made nonsingular and chosen 
to be the same as that for the field strength^ It then follows that the 
corresponding contribution D- to the index D is vanishing. 

In consequence, the index D of (|C.ip is equal to -|-1. Immediately, one 
notices that this index D is nothing but the index of the singular point of the 
wavefunction g'_(k). The importance of such singularities of wavefunctions 


^Restricting the region to one-half of the Brillouin zone is essential to this argument, 
because one cannot hnd a nonsingular connection on the whole Brillouin zone. 
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on the Brillouin zone has already been pointed out by Kohmoto in his early 
work m . Furthermore, one can easily see that this index D is equal to the 
present Z2 index of (13.41) by recalling the argument of Appendix iBl 
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